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Abstract 

^ I Long-range contributions to the three-nucleon force that have been recently worked out in chiral 

p^ I effective field theory at next-to-next-to-next-to-leading order are for the first time included in the 

triton and the doublet nucleon-deuteron scattering length calculations. The strengths of the two 
short-range terms available at this order in the chiral expansion are determined from the triton 
binding energy and the neutron-deuteron doublet scattering length. The structure of the resulting 
\^ ', three-nucleon force is explored and effects for the two-nucleon correlation function in the triton are 

investigated. Expectation values of the individual contributions to the three-nucleon force in the 
triton are found to be in the range from a few 100 keV to about 1 MeV. Our study demonstrates that 



> 



in 



(~^ • the very complicated operator structure of the novel chiral three-nucleon forces can be successively 



% 



implemented in three-nucleon Faddeev calculations. 
PACS numbers: 21.45.-v, 21.30.Fe, 21.45.Ff 



I. INTRODUCTION 

Chiral effective field theory (EFT) provides a powerful framework to systematically de- 
scribe low-energy dynamics of few- and many-nucleon systems. Various variants of effective 
theories for nuclear forces have been explored, see [l|-l3| for recent review articles. Up to now, 
the most advanced few-nucleon studies have been carried out within a framework based on 
pions and nucleons as the only explicit degrees of freedom taken into account in the effective 
Lagrangian. Within this approach, the nucleon-nucleon (NN) force is currently available 
up to next-to-next-to-next-to leading order (N^LO ) in the chiral expansion. At this chiral 
order, it receives contributions from one-, two- and three-pion exchange diagrams as well 
as short-range NN contact interactions with up to four derivatives. As demonstrated in 
Refs. yi |5[, NN phase shifts are accurately described at N^LO up to laboratory energies of 
the order of 200 MeV. The theoretical uncertainty due to truncation of the chiral expan- 
sion is estimated in Ref. j^ by means of a cutoff variation. Within the spectral function 
regularization (SFR) framework [6] adopted in Ref. jj], the NN potential depends on two 
ultraviolet cutoffs A and A. The first one removes large-mass components in the spectrum of 
the two-pion exchange potential which cannot be correctly described within the chiral EFT 
framework while the other one provides regularization of the Lippmann-Schwinger equation. 
Five different combinations of these cut-off parameters are available for the NN potentials 
of Ref. [^ . The residual dependence of low-energy observables on the cutoff choice provides 
a measure of the importance of higher-order contact interactions and thus may serve as an 
estimate of the theoretical uncertainty. 

Parallel to these developments three-nucleon force (3NF) has also been explored within 
the framework of chiral effective field theory. The first nonvanishing contributions to the 
3NF emerge at next-to-next-to-leading order (N^LO ) p] from the two-pion exchange and 
one-pion-exchange-contact diagrams as well as the purely short-range derivative-less three- 
nucleon contact interaction [TJ, see also Ref. J8| for a pioneering work along this line. The 
resulting N^LO three-nucleon potential depends on two low-energy constants (LECs) D and 
E accompanying the short-range nNN and NNN vertices, respectively. The values of these 
LECs need to be fixed from a fit to few-nucleon data. Among a few possible observables that 
have been used in this connection are the triton binding energy and the nucleon-deuteron 
doublet scattering length ^a„rf |7|,l9|, a-particle binding energy jloL 11|. the properties of light 



nuclei [l2| and the triton beta deca y |l3| . The N^LO 3NF of [7| was successfully used in 
three-body calculations, see Refs. [ij, |l5[ for a few examples of recent studies. At this order, 
chiral EFT yields a good description of elastic scattering and deuteron breakup observables 
up to energies of about ^ 50 MeV. The accuracy of the results in this regime is comparable 
with the one that is achieved by realistic phenomenological NN and 3N interactions such 
as e.g. AV18 fiel 2NF in combination with UrbanalX Jl^LSNF or CD-Bonn fisl 2NF in 
combination with the Tucson- Melbourne [19|] 3NF, see [14, |20|. However, the spread of the 
results is relatively large for some spin observables which clearly calls for the inclusion of 
new terms of the nuclear interaction that occur at higher orders of the chiral expansion. 

Subleading contributions to the 3NF are currently being investigated by several groups. 
At N^LO , one has to take into account (irreducible) contributions emerging from all possible 
one-loop three-nucleon diagrams constructed with the lowest order vertices. In addition. 



there are (tree-level) leading relativistic corrections, see [2l| for an early work on the longest- 
range relativistic corrections. Note that the tree diagrams involving higher-order vertices 
from the effective chiral Lagrangian do not produce any irreducible pieces. Effects due 



to two-pion exchange 3NF in elastic nucleon-deuteron scattering were already explored by 



Isliikawa and Robilotta [22| within a hybrid approach and found to be rather small. The 
N'^LO contributions feed into five different topologies as will be explained in detail in the next 
section. The explicit expressions both in momentum and in coordinate space for the long- 
range contributions have already been worked out 22, l23| . Their inclusion in numerical few- 
body calculations appears to be challenging due to the very rich and complicated operator 
structure. The large number of terms in the 3NF at N^LO , see Ref. |23j, requires an efficient 
method of performing the partial wave decomposition. Recently such a method has been 
proposed [2J] and tested for the Tucson-Melbourne force [25| . Here and in what follows, this 
approach will be referred to as the automatized partial wave decomposition (aPWD). In this 
paper we apply this method of the numerical partial wave decomposition to the N^LO 3NF 



contributions derived in 23|]. For the first time, the parts of 3NF at N'^LO different from 
the two-pion exchange force are included in the triton and the scattering length calculations. 
In order to test the implementations and get a first hint to possible effects of these forces, 
we fix the two LECs entering the 3NF from the triton binding energy and the nucleon- 
deuteron doublet scattering length and explore the effects due to these novel 3NF terms by 
computing the ^H properties. Although this calculations is still incomplete since not all 3NF 
contributions at N^LO are taken into account, it provides an important first step towards 
the complete N'^LO analysis of 3N scattering and demonstrates our ability to numerically 
handle the rather complicated structure of the subleading chiral 3NF. 

Our paper is organized as follows. In Sect. [ITl we describe briefiy the structure of the 
chiral 3NF at N^LO . In Sect. Ulllwe discuss in detail the partial wave decomposition needed 
in our scattering and bound-state calculations. Next, the procedure of fixing the LECs is 
described in Sect. HVl where also the obtained values of LECs are listed. These results are 
used in Sect. |V]to explore some properties of the triton. Finally, our findings are summarized 
in Sec. |VT1 

II. 3NF AT N^LO 

The subleading (i.e. N^LO ) contributions to the three-nucleon force V123 can be written 
in the form 23 



V^123 — V21T + V2tt-1tt + Vring + V^lTT-cont + V2TT-cont + V^l/m • (2-1) 

where the individual terms refer, in order, to the two-pion exchange, two-pion-one-pion- 
exchange, ring (i.e. one pion being exchanged between each of the three nucleon pairs), 
one-pion-exchange-contact and two-pion-exchange-contact contributions as well as the lead- 
ing relativistic corrections, see Fig. 1 of Ref. 23j] for a diagrammatic representation. The 



expressions for the (static) long-range part of the 3NF given by the first three terms in the 



above equation have been worked out in heavy-baryon chiral perturbation theory in Ref. |23 . 
The two-pion exchange contribution at the one-loop level has also been calculated within 
the infrared-regularized version of chiral perturbation theory in Ref. [22| . The shorter-range 
contributions involving two-nucleon contact interactions and relativistic corrections are cur- 
rently being worked out [261] . 

While the two-pion-exchange V2n and one-pion-exchange-contact parts ViT,-cont already 
occur at N^LO and receive corrections at N^LO , the remaining topologies first emerge 
at N'^LO . It is important to emphasize that all subleading contributions to the 3NF are 



parameter- free. Thus, the low-energy constants D and E entering the one-pion-exchange- 
contact and the purely short-range parts of the 3NF at N^LO are the only unknown param- 
eters up to N^LO . 

Here and in what follows, we adopt the notation in which a given 3NF \/i23 is decomposed 
into three terms 

V,,, = \/(^) + \/(2) + y(3) ^ (2.2) 

where each \^^*^ is symmetrical under interchanging the nucleons j and k {i,j,k = 1,2,3, 
i ^ j ^ k). Clearly, this condition does not specify V'^^^ uniquely. In the following we choose 
V^^^ in such a way that the number of operator structures is minimized which is convenient 
for the aPWD. 

The operator structure of the 27r exchange part V2J at N^LO remains the same as at 
N^LO 

r£^ = Fi Cf2 ■ q2 0^3 ■ gs T2 ■ T3 + F2 0^2 ■ ^2 0^3 ■ ^"3 02 >^ Qs ' (^l T2 X Tg ■ Ti, (2.3) 

where ql is the momentum transfer to the i-th nucleon, ql + q2 + Qs = and (fi (Ti) are 
Pauli spin (isospin) matrices for nucleon i. The scalar functions Fi = Fi{q2,q3,q2 ■ qs) and 
F2 = F2{q2,q3,q2 ■ qs) depend on the LECs £1,3,4 which accompany the subleading pion- 
nucleon vertices. Chiral expansion of Fi and F2 up to N^LO has the form |23| 

p ^ A (-4ciM,^ + 2d3g-2 ■ gi) , p ._. 

' 4F4 (gi + M2)(g| + M2) + ^ ^- > 

- - ^^ ""' +F2, (2.5) 



4F^ {ql + M^^){ql + M^^) 

with 

g\ 1 

1287rF6(gi + M2)(g| + M2) 

+ {2Ml + ql + ql + 2g-2 ■ (b){^Ml + 3^2' + H + 4gi ■ (b)A{\q2 + gi|)) (2.6) 



^1 = .oo_^./ 2 , ,^2^/ 2 , ...^ (^^(^' + 3g2' + 3gI + 4gi•gi) 
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(M^ + (4M^^ + qi + gg' + 2g-2 ■ q:,)A{\q2 + gi|)) , (2.7) 



1287rF6(g| + M2)(g| + M2) 
where the loop function A{q) is defined as 

Aiq) = — arctan— ^. (2.8) 

^^^ 2g 2M^ ^ ^ 

The axial- vector coupling constant, the weak pion decay constant and the pion mass are 
denoted as qa, F^, and M^, respectively. Note that SFR changes A{q) (see |6|). For the 
study here we do not need to consider this change since they differ only by higher order 
polynomials [l]. The quantities q appearing in the above expressions are related to the 
N^LO LECs Cj entering the effective chiral Lagrangian via 



where we have adopted the values for the q from Ref. jj|, namely 

ci = -0.81 GeV-\ C3 = -3.40 GeV"\ C4 = 3.40 GeV"^ . (2.9) 



These values are consistent with the ones determined from pion-nucleon scattering [27| . The 
finite shifts of the LECs q in Eq. f l2.9p emerge from pion loops at N^LO . 
The V^TT-ivr interaction at N^LO has the following operator structure 



T/(l) _ ^1-gl 

Vn 



T2 ■ Ti (^3 ■ g2 92 ■ qi Fi{q2) + <T3 " ^2 ^^2(^2) + ^Z ' Ql ^^3(^2)) 



2^-1" ql + M^V 

+ T3 ■ Ti (ct2 ■ g2 92 ■ qi -P4(g2) + <?2 ■ gi -F5('?2) + 0^3 ■ g2 -F6('?2) + 0^3 • ?! -^7(52)) 

+ T2 XT3-ri(T2 Xa3-g2i^8(g2) +(2^3). (2.10) 

with scalar functions Fi{qi). The interchange of nucleons (2 o 3) refers to the interchange 
of momentum vectors, spin and isospin matrices and arguments of the functions Fj. Ex- 
plicit expressions for the scalar functions Fi{qi) which appear in the Vg^-i^ can be found in 
Ref. H. 

The Vjjg force is chosen as 

King = O^i ■ 0^2 T2 ■ T3 i?i + (Ti ■ ^1^2 " ?! T2 " T3 i?2 + <?! " giO'2 " ^3 ^2 " T3 i?3 

+ <?1 ■ 'f3'?2 ■ gi T2 ■ T3 /?4 + <?1 ■ g3'?2 " 53 T2 " T3 R5 + Ti ■ T3 i?6 + <?! " 'fl'?3 " Ql Rj 
+ (Ti ■ gi(T3 ■ g3 i?8 + 0^1 ■ g3'?3 " gl i^g + O^l " '?3 ^10 + gi ■ gs X 0^2 Ti ■ T2 X T3 i?ii 

+ Ti-T2Si + (Ti- qia-i • gl Ti ■ r2 5*2 + CTi ■ g3o'3 • gi Ti ■ T2 53 + ai ■ gia3 ■ q-^Ti- T2 ^4 

+ <?! ■ g3'?3 ■ gs Ti ■ T2 55 + ai ■ CT3 Ti " T2 5*6 + gl ■ g3 X ^i Ti ■ r2 X T3 ^t 

+ (2^3), (2.11) 

where the expressions for the scalar functions Ri = -Rj(gi, g3, qi ■ qs) and Si = Si{qi, q^, gi ■ g'3) 
can be found in Ref. 23 . 

The modifications of the one-pion-exchange-contact V^iV_cont term arising at N^LO are in 

preparation [26]. Thus, instead of the full V^i^_cont interaction we use the lowest-order result 
for it resulting at N^LO fl\: 

X .(1) 9a D ai ■ qi - - , - - \ tn ^n\ 

Vd-term = 'J^ ^2 ^ j^2 ("^1 " ^3 CT3 " gi + Ti ■ T2 ^2 " gi) . (2.12) 

The low energy constant D can be expressed as D = c/)/(i^^A^), where c^ is a dimensionless 
free parameter and the chiral symmetry breaking scale A,^ is estimated to be A^ = 700 MeV. 
Here and in what follows, we use F^ = 92.4 MeV for the pion decay constant. The value 
of cd has to be determined from experimental data. This is described in section |IV] for our 
fit for the test case of this study. Of course, results of this fit will significantly change when 
the complete short range interaction is taken into account. 

As already pointed out above, the remaining terms V2j_cont ^^^1 V^l^ are also not available 
yet |26j and thus cannot be taken into account in the present study. Finally, the purely 
short-range part of the 3NF has the form [2] : 

V^Xrm = ET2-T,. (2.13) 



Again, the LEG E is usually expressed in terms of a dimension-less parameter ce via E = 
ce/{E^A^) which needs to be determined from (at least) three-nucleon data. 

To summarize, in this paper we use the V2J, y-in-i-K^ ^ring N'^LO terms combined with 
the V^ilL.™ and V^Xrm terms at N^LO. The N^LO contributions to viZ^^, and Vi^_^^^^ 
and relativistic corrections V^^^ are not included. Finally, the remaining terms l^*-^^ and 
V'^^'^ of Eq. (12. 2p can be obtained from V^^'^ by appropriate permutations of nucleons. 



III. NUMERICAL CALCULATIONS OF 3NF MATRIX ELEMENTS 

We work in the momentum space using three-nucleon partial-wave states | p, g, a) in the 
jj-coupling [281. 1291 



I p,g,«) =1 vq{ls)3{^\)lUl)JMj) I (ti)rMT) , (3.1) 

where p and q are magnitudes of the standard Jacobi momenta and a denotes a set of discrete 
quantum numbers defined in the following way: the spin s of the subsystem composed from 
nucleons 2 and 3 is coupled with their orbital angular momentum / to the total angular 
momentum j. The spin 1/2 of the spectator particle 1 couples with its relative orbital 
angular momentum A to the total angular momentum / of nucleon 1. Finally, j and / are 
coupled to the total 3N angular momentum J with the projection Mj. For the isospin part, 
the total isospin t of the subsystem (23) is coupled with the isospin 1/2 of the spectator 
nucleon to the total 3N isospin T with the projection Mt- 

The matrix elements of V2J ■, V27r_i7r and V^^J^^ forces in the basis | p, g, a) are obtained us- 



ing the recently proposed aPWD method |24,l25|]. In this approach the spin- momentum and 



isospin parts of three-nucleon interactions are calculated using a software for symbolic cal- 
culations. The resulting momentum-dependent functions are then integrated numerically in 
five dimensions over angular variables. The major advantage of this method is its generality. 
It can be applied to any momentum-spin-isospin operator including, in particular, the full 
operator structure of the 3NF at N^LO including even the non-local relativisitic corrections. 
The only complication emerges in the treatment of the ring contributions to the 3NF due to 
rather complex expressions for the functions Ri = Ri{qi, qs, qi ■ qs) and Si = Si{qi, qs, qi ■ gs) 
which involve certain scalar integrals related to the three-point function. For any given val- 
ues of the arguments, these integrals have to be computed numerically. This evaluation is 
too expensive to be carried out on-the-fly during the aPWD. Moreover, while the functions 
Ri and Si are, of course, finite and smooth for all possible values of their arguments, they 



are given in Ref. |23[ as linear combinations of terms, some of which are becoming singular 
under certain kinematical conditions. Their numerical implementation therefore requires 
special care. In order to deal with these difficulties, we first evaluate the functions Ri and 
Si at a (fixed) dense grid of points for their arguments and then use interpolation to com- 
pute them for arbitrary values of gi, gs and qi ■ gs as needed in the aPWD approach. We 



have carefully checked the stability of this procedure by increasing the density of the grid 
points. Finally, the partial wave decomposition of VJ-_jg^^ and V}_{^^^ is performed with the 
standard techniques |7| but also verified with the new method. 

Examples of the resulting matrix elements {p', q', (y'\V,- \p, g, a) are given in Fig. [T] as a 
function of the momentum p. Here, we fix the momenta to be p' = 0.268 fm~^, g' = 2.842 
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TABLE I: The values of the discrete quantum numbers for the selected a states (13. ip with the 
total angular momentum J = 1/2, the total isospin T = 1/2 and its projection Mt = —1/2. 



fni~^ and q = 0.132 fni~^ and consider the following four channel combinations: {a' = a = 
1), {a' = l,a = 5), {a' = a = 4) and (a' = 4, a = 7). These channels correspond to the 
quantum numbers given in Tab. [H Channels a = 1 and a = 5 are especially important 
since those two states provide the dominant components of the ^H wave function. As can 
be noticed, all three components of the N^LO 3NF give a strong contribution: the VjJ„ 

1, a = 5) and V^^ 



1), ^£-- for ^«' 



dominates for the channel combination {a' = a 
for {a' = a = 4) and {a' = 4, a = 7). We emphasize, however, that the large size of these 
matrix elements (which contain certain admixtures of short-range operators) as compared 
to the N^LO terms does not necessarily imply that theirs effects in low-energy observables 
are large. The values of the LECs for V^_[^^^ and V}_{^^^ terms shown in Fig. [T] are set 
to he cn = 1 and c^; = 1 in order to allow for a qualitative comparison of the strength of 
the individual terms. Their real contributions emerging after fitting the LECs cd and ce to 
experimental data will be discussed in the next section. Fig. [1] also clearly demonstrates that 
not all terms contribute to each channel combination due to the spin-isospin dependence, 
o,^^ T/v terms contribute to V^-^^ only to {a' = l,a = 5) and the V^ 



2n-lTT 



The V,^:U^ and V,_,,_ 
term only contributes to (a' 



(1) 

e—term 



a = 1). Further, V^jJ^ vanishes for a' 



a 



4. 



The cornponent V^^^ of V^i23 enters the dynamical equations for 3N bound and scattering 
states 29|, l30|] (see below) only in combination with the permutation operator P forming 
the operator V^^\l + P). The permutation operator P = P12-P23 + -P13-P23 is built from the 
transpositions Pij, which interchange nucleons i and j. The aPWD scheme can be used to 
obtain directly the ^'•^^(1 -|- P) matrix elements 25| which allows us to avoid uncertainties 
associated with the partial wave decomposition of the permutation operator. The resulting 
matrix elements {p',q',a'\V^^\l + P)\p,q,a) are given in Fig. [2] for the same momenta 
and channel combinations as in Fig. [T] Again, we observe that matrix elements of all 
contributions to the 3NF are large for these momenta. The Vj„g contribution is particularly 



(1) 



Ki) 



large in matrix elements with a' = 1. For these channels also V2J and Vg^r-ivr ^^^ significant. 
Moreover for a' = a 
and Ce = 1- For a' = 



Ki) 



Ki) 



1 also V^4^^.^ and V^'_t^^^ 



are non- negligible, at least with c/5 = 1 
4, the V2J piece dominates. Nevertheless in the case of a = 7 also 



^ring P^ovides a significant contribution. Due to their spin-isospin structure the V^ L and 

K_L„^ forces are absent for (a' 



^(1) 

ring P^^vnaco a oij^iiiin.aiit ^.wiitiiuutiwii. j^ uc tw uiicii opiii-iowopiii otiui^tuic uiic v d-term 

/^e-term fo^ccs are abseut for {a' = a = 4) and {a' = 4, a = 7) channel combinations. 

The 27r-exchange contribution at N^LO can be compared with a corresponding part of the 
N^LO interaction. The example is given in Fig. [3l where the dashed (solid) line represents 
the predictions obtained at N^LO (N^LO ). The Vg^^ (1 + P) niatrix elements differ signif- 
icantly for all considered channel pairs. Fig. [3] demonstrates that these differences mainly 
originate from the parts of the V2J interaction proportional to the Fi and F2 formfactors, 
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1: (Color online) The N^LO 3NF matrix elements (before regularization) 
'\V^^'\p,q,a > as a function of the momentum p for p' = 0.268 hji~^,q' = 2.842 
fm^^ and q = 0.132 fm~^ and different combinations {a', a): a) (1,1), b) (1,5), c) (4,4) and d) 
(4,7). The solid (black), dashed(red), double-dash-dotted(green), dotted(blue) and dash-dotted 
(magenta) line shows the V27r,V27r-i7r,K-mg,Vrf-term and Ve-term components, respectively. 



see Eqs. (l2.4p - (I2.5I) . In particular, in Fig. |3l we also show results for V^^ (l + P) at N'^LO but 
with Fi and F2 being artificially set to zero. In such a case the only difference between 
matrix elements at N^LO and N^LO comes from the different values of q and q LECs. 
The presented matrix elements have then a similar dependence on momentum p as well as 
a similar magnitude in N^LO and N^LO . Notice, however, that these observations do not 
necessarily imply that the N^LO corrections to the 3NF lead to large effects in low-energy 
three- nucleon observables. In fact, the opposite was observed in Ref. 22| for the case of the 
two-pion exchange topology. 

The V^ and V^ (1 + P) matrix elements shown in Figs. [1] and [2] have to be regularized 
prior to being used as input to the dynamical equations [Tj . We use the regulator function 
of the form [l| 

f{p, q) = exp — — (3.2) 



(4A2)3 

which ensures that the large momenta are sufficiently suppressed. Following Ref. (l| we use 
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FIG. 2: (Color online) The N'^LO 3NF matrix elements (before regularization) 
< p',q',a'\V^^'{l + P)\p,q,a > as a function of the momentum p. The momenta p',q' and 
q, channel combinations (q',q) and lines are the same as in Fig. [TJ 



three values of A parameter: 450, 550 and 600 MeV. The regularization transforms matrix 
elements as 



{p',q',a'\V^'\l + P)\p,q,a) ^ f{p',q'){p',q',a'\V^'\l + P)\p,q,a)f{p,q). 



(3.3) 



The examples of the regularized V^^\l + P) matrix elements are compared to the nonreg- 



ularized ones in Fig. HI This is done separately for V2V (1 + -P) and V)^'{l + P) contributions. 



(1)/ 



^(1) 



The momenta are p' = q' = q = 0.132 fm~^ (upper row) and p' = 0.268 fm~^, q' = 2.842 
fm~^ and q = 0.132 fm~^ (lower row) and we show only the {a' = a = 1) channel combina- 
tion. In the upper row, where the momenta p', q' and q are small, all regulator functions are 
close to 1 for small values of p. For momenta p > 1 fm~^, the different A values lead to dif- 
ferent slopes of matrix elements. The lowest value of the parameter A = 450 MeV forces the 
fastest decreasing of V-^ (1 + P) matrix elements. In such a case, the short-range part of the 
interaction is suppressed. On the contrary, the highest value A = 600 MeV allows for larger 
contributions of short-range interactions. In the lower row, where the momenta p' and q' are 
bigger, the effects of the regularization are seen already at the low values of p. For p = 0.001 
fm~^ the regularization factor f{p', q')f{p, q) changes from 0.194 for A = 450 MeV to 0.747 
for A = 600 MeV. This strong cutoff dependence is expected to be largely compensated by 
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FIG. 3: (Color online) The 27r-exchange part of 3NF matrix elements (before regularization) 
< p', q', a'l V^^ (1 + P) \p, q,a > as a fmiction of the p momentum. The momenta p', q', q and chan- 
nel combinations (a', a) are the same as in Fig.[TJ The dashed (black) and solid (red) lines represent 
N^LO and N^LO results, respectively. The dotted (black) line describes the matrix elements of 
V^^ (1 + P) obtained with Fi and F2 set artificially to zero (see text). 

an appropriate "running" of the LECs cd and ce when calculating low-energy observables. 



IV. DETERMINATION OF THE LECS cd AND ce AT N^LO . 

Once the new terms are added to the 3NF, the procedure of refitting of the LECs cd and 
Ce has to be repeated. We follow Ref.[7| and use the triton binding energy E ^ and the 
neutron-deuteron doublet scattering length "^and as two observables from which cd and ce 
can be obtained. The up-to-date experimental values are E ^ = 8.481821(5) MeV [31] and 



'^and = 0.645(7) fm [32 



Our procedure to fix the values of LECs can be divided into two steps. First, the depen- 
dence of £" ^ on Ce for a given value of c^ is determined. The requirement to reproduce 
the experimental value of the triton binding energy yields a set of combinations cd and ce- 
This set is then used in the calculations of ^a^^ what allows us to find which pair of cd and 
Ce describes both observables simultaneously. Such a procedure has to be repeated for all 
A-values used in the regularization. The same values of the cutoff A are used to suppress 
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FIG. 4: (Color on-line) The N^LO 3NF matrix elements < p',q',a' = 11^/^1 + ^)b,9,a = 1 > 
as a function of the p momentum for the momenta p' = q' = q = 0.132 fm^^ (top) and p' = 0.268 
im-^,q' = 2.842 fm'^ and q = 0.132 fm"! (bottom). The two components of the N^LO 3NF are 



.(1), 



.(1) 



(1)/ 



shown: V^!,' {I + P) (left) and T/;4g(l + P) (right). The sohd (black) line represents 1^2^^ (1 + P) 

.(1) 



and Vmg(l + ^) niatrix elements before regularization. The dashed(red), dash-dotted (green) and 

d(bh 
in dSSD with A = 450, 550 and 600 MeV, respectively. 



dotted(blue) curve represents the ^2,^ (1 + P) and Vmg(''^ ~^ ^^ matrix elements regularized as 



high momenta in the NN potential in order to ensure the convergence of the integral in the 
Lippmann-Schwinger equation. The chiral NN potential depends, in addition, on another 
cut-off parameter A emerging from the SFR of the two-pion exchange potential. We follow 
Ref. [ll and use five combinations of (A, A) shown in Tab. HIl 

We compute the '^H wave function using the method described in [30]- Here we mention 
only that the full triton wave function ^ = {1 + P)tp is given by its Faddeev component i/) 
being the solution of the Faddeev equation 



^ = GotPi; + (1 + Got)GoV^'\l + P)^ . 



(4.1) 



Here, Gq is the free 3N propagator, P is the same permutation operator as defined above and 
t is the two-body t-operator generated from a given NN potential through the Lippmann- 
Schwinger equation. 

We use the 3N states |p, g, a) defined on the grids of 68 p-points and 48 g-points in 
intervals p G (0, 15) fm~^ and q G (0, 10) fm~^, respectively. We take into account all states 



11 



3.00 



"1 — I — I — I — I — I — I — I — I — I — I — I — I — r 




0.00 I — ^ I I I I I I I I I I I I L^ 




D 



FIG. 5: (Color online) The intermediate results obtained during the fitting procedure for the cd 
and ce N^LO LECs for the 3rd cut-off combination from Tab. [Til the values of the cd and ce 
LECs which give the experimental E ^ (top) and the dependence of ^a„rf on cd (bottom). The 



dotted line shows the experimental value ^a„d=0. 645(7) fm [32 1. 



up to the two-body total angular momentum j = 5 for the NN potential, and all states up 
to j = 3 for the 3N interaction. 

We solve Eq. (14. ip and find pairs of the LECs cd and ce which reproduce the experimental 
value oi E ^ . It is exemplified in the top panel of Fig. |5] for the 3rd cut-off combination 
from Tab. [TTl The dependence is smooth and for some values of ce there are two possible 
values of cjj. 

In the second step of the fitting procedure the doublet scattering length ^a^^ is calculated 
for the (cd, ce) pairs, which reproduce the correct value oi E ^ . To this aim we first solve 
the Faddeev equation for the auxiliary amplitude T at zero incoming nucleon energy, 

T = tP(P + (l + tGo)V«(l + P)0 + tPGoT + (l + tGo)V«(l + P)GoT, (4.2) 
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cut-off 


(A, A) 


CD 


ce 


1 


(450,500) 


10.78 


-0.172 


2 


(600,500) 


12.00 


1.254 


3 


(550,600) 


11.67 


2.120 


4 


(450,700) 


7.21 


-0.748 


5 


(600,700) 


14.07 


1.704 



TABLE II: The values of cd and ce LECs for the different parametrizations of chiral N^LO po- 
tential. 

where the initial channel state occurring in the driving terms is composed of the deuteron 
and a plane wave state of the projectile nucleon. The amplitude for the elastic nucleon- 
deuteron scattering is then given by 



U 



PGo' 



PT + 1/(^)(1 + P)0 + \/(^)(i + p)Gor. 



(4.3) 



We refer to [29|, |33| for a general overview on 3N scattering and for more details on the 
practical implementation of the Faddeev equations. The expression for ^a„rf in our basis and 
further technical details can be found in [3J]. In this second step of the fitting procedure we 
use grids of 32 p points in the range p G (0, 25) fm~^ and 31 q points in the range q G (0, 15) 
fm~^. Similarly to the triton calculations, the NN (3N) potential acts in all states up to 
j = 5(3). Our calculations are accurate up to 2 keV for the binding energy and up to 
0.005 fm for the scattering length. We expect that, even in N^LO , the chiral expansion 
of the nuclear forces induces uncertainties that are larger than these estimates. Therefore, 
the numerical calculations are sufficiently accurate to perform sensible fits of cd and ce at 
N^LO . 

The final values of cd and ce LECs, which reproduce the experimental values of E ^ and 
'^ttnd are given in Tab. [TTl For all combinations of cut-off parameters, the LEG cd remains 
positive with the value around 10. It weakly depends on the value of the cutoff A and 
becomes larger with increasing A. The second LEG ce changes in a more complicated way. 
While again it is smallest in magnitude for the smallest value of A, its biggest value is for 

the medium A=550 MeV and then decreases while moving to A=600 MeV. Note, that ce 

(1) 



changes sign so the V}_1^^.^ interaction changes from attractive to repulsive. We also stress 
that while the value of the LEG c/? appears to be rather large, the expectation value of 
the one-pion-exchange-contact part of the 3NF is of a natural size. It remains to be seen 
whether the complete calculation including the remainig 3NF contributions at N^LO will 
lead to more natural values of the LEG cd- 

Finally, we would like to emphasize that the values of the LEGs which are bare pa- 
rameters must be refitted at each order in the chiral expansion (and, of course, for each 
cutoff combination). This is in contrast with chiral perturbation theory calculations in the 
Goldstone-boson and single-nucleon sectors where the scattering amplitude is usually ex- 
pressed in terms of renormalized LEGs. Indeed, using the values of cn and ce determined at 
N^LO in the N^LO calculation would generally result in a poor description of low-energy ob- 
servables. For example, for the 3rd cutoff combination, the 3NF at N^LO furnished with the 
N^LO values cd = -0.45 and ce = -0.798 [l| yields E'" = -8.197 MeV and ^and = 1-004 
fm which are far from the experimental values. Similarly, while the combination cd = 1.5744 
and Ce = —17.8 allows to reproduce the triton binding energy and the doublet nd scattering 
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cut-off 


{Ho) [MeV] 


{Vnn) [MeV] 


{Vsn) [MeV] 


1 


35.972 


-43.459 


-0.994 


2 


54.708 


-61.515 


-1.673 


3 


48.088 


-55.187 


-1.381 


4 


33.232 


-41.050 


-0.663 


5 


53.504 


-60.278 


-1.706 



TABLE III: Expectation values {Hq), {Vnm) and (Vsn) in the triton for different parametrizations 
of tfie chiral N^LO potential as discussed in the text. 

length for N^LO NN force accompanied with the N^LO 3NF, it produces E^^ = -7.542 MeV 
and "^and = 1.4354 fm when the N^LO 3NF is used. Note that a big value of ce obtained in 
this case seems to violate naturalness. These results demonstrate clearly that fitting to the 
data has to be made consistently within the given order of the chiral expansion. Therefore, 
also the results of our test fit here, which does not include the full 3NF at N^LO , has to be 
taken with care. However, the fit results allows us to study the properties of ^H in the next 
section. 



V. THE PROPERTIES OF =^H 

Once the values of cn and ce low energy parameters are established, one can explore the 
properties of the ^H wave function. We begin with the expectation values of the kinetic 
energy {Hq), the NN potential energy (V^n) and the 3N potential energy {Vsn) which are 
listed in Tab. IIIII The expectation values clearly depend on the cut-off parameters A and 
A as they should. Not surprisingly, the expectation values of both the NN potential and 
the 3NF are smallest for the softest cutoff A = 450 MeV. Higher A-values lead to stronger 
3NF contributions to the E ^ which for A = 600 MeV reaches about ~ 20%. For all cutoff 
values, one clearly observes the dominance of the NN forces, (V3Ar)/(V7VAr) = 1.5... 3%, 
in agreement with the expectations based on the chiral power counting. Note that the 
interplay of N^LO counter terms and N^LO structures of the 3NF does not allow to use 
these expectation values for an assessment of the contributions of N^LO 3NFs. Similarly 
as in the NN interaction, one observes strong cancelations of N^LO contributions also with 
N^LO contact interactions. Such a comparison only makes sense for renormalized quanties, 
which we are not able to identify here. One also observes that the dependence of the 
expectation values on the SFR cutoff A is less pronounced as the A-dependence. For example 
for cutoff combinations 1 and 4, which differ only in the choice of A, the (Hq) and (Vnn) 
differ by about 3 MeV. On the other hand, the differences reach almost 20 MeV for the 
cutoff combinations 1 and 2 which have the same SFR cutoff A but different values of A. 
This holds true also for (V^n). 

A more detailed information about the 3NF triton expectation values is given in Tab. IIVI 
The expectation value of the 3NF (Vsat) is split into the individual contributions from various 
topologies. The expectation value of the two-pion-exchange potential {V2n) shows a smooth 
dependence with A. For the softest cutoff A = 450 MeV, the two-pion exchange 3NF 
turns out to be most attractive providing more than 0.5 MeV to the triton binding energy. 
With increasing A, the contribution of V27r becomes weaker. For the cutoff combination 2, 
A = 600 MeV, A = 500 MeV, the additional binding due to the two-pion exchange 3NF 
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cut-off 


(^2.) [MeV] 


(^2.-1.) [MeV] 


{Vring) [MeV] 


{Vd-term) [MeV] 


{V,_term) [MeV] 


1 


-0.639 


0.458 


-0.147 


-0.693 


0.027 


2 


-0.241 


-0.580 


-1.114 


0.694 


-0.432 


3 


-0.473 


0.107 


-0.191 


-0.708 


-0.116 


4 


-0.771 


0.539 


-0.452 


-0.259 


0.281 


5 


-0.377 


-0.275 


-0.622 


-0.119 


-0.313 



TABLE IV: The expectation values for the different parts of the 3N potential and for the different 
parametrizations of the chiral N^LO potential. 

only amounts to about 240 keV. Interestingly, most of the attraction necessary to reproduce 
the triton binding energy is produced in this case by the ring topology, which is found to 
be attractive for all cutoff combinations. Contrary to the longest-range two-pion exchange 
topology, the contributions of the ring diagrams are enhanced for the largest value of the 
cutoff A = 600 MeV. Qualitatively, this behavior might be expected given the fact that 
the large values of A probe the shorter-range part of Vrmg which is of the van der-Waals 
type, i.e. the matrix elements grow rapidly with decreasing relative distances between the 
nucleons. The cutoff dependence of this contribution is very large. This explicitly shows the 
dependence on the short distance pieces making it impossible to estimate the impact of this 
topology on low energy observables based on our results. The expectation value of the two- 
pion-one-pion-exchange topology (y2-K~i-K) also strongly depends on A. It changes sign from 
positive for A = 450 MeV to negative at A = 600 MeV. The dependence on the SFR cutoff 
is fairly weak. Note however that a stronger dependence might be induced once the SFR 
regularized A{q) has been taken into account. The (yd-term) shows the most complicated 
behaviour. It achieves the lowest value for the intermediate value of A (cutoff combination 
3) and, for A = 600 MeV, shows a strong dependence on the SFR cutoff A. In particular, 
for the smaller value A = 500 MeV, this contribution to the 3NF becomes repulsive and 
relatively big, while for A = 700 MeV the expectation value remains negative. Finally, the 
(ye-term) expectation value changes smoothly with A. It also changes its sign from positive 
at the lowest A to negative at A = 600 MeV. Last but not least, we emphasize that the 
expectation values discussed above as well as the separation of the potential energy into the 
contributions due to the NN potential and 3NF do not correspond to observable quantities 
and are expected to show a strong cutoff dependence. Notice further that expectation values 
of the various 3NF contributions are, strictly speaking, bare quantities. Comparing their 
size with the one of the N^LO terms does, therefore, not allow to draw conclusions about 
the convergence of the chiral expansion [38|. It is comforting to see that all expectation 
values turn out to be of a reasonable size. 



We now turn to the two-nucleon correlation function of '^H which is defined as 30| 



C(r) 






(5.1) 



Here, r7j is the relative distance corresponding to the Jacobi momentum p. In Fig. |6] the 
correlation function is shown for the same combinations of the regularization parameters 
as in TabJTTl Thin lines represent predictions based on NN interactions only while thick 
ones show the predictions based on NN-I-3N forces. As expected, the softest cutoff value 
A = 450 MeV yields a more flat correlation function with the less amount of the short-range 
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FIG. 6: (Color online) The two-body correlation function for the triton for different (A, A) pairs 
given in Tab. HIl The thin dotted (black), dashed (red), solid (blue), dash-double-dotted (black) 
and dash-dotted (red) lines correspond to predictions based only on NN interaction with cut-off 
numbers from 1 to 5, respectively. The thick lines represent predictions with the same regularization 
parameters but based on NN and 3N forces. 

correlations. The higher A values prefer distributions concentrated around the maximum at 
r f« 1.5 fm. The effects of the 3NFs are small for the lowest A but increase with increasing 
A. For the lowest A there is also a strong dependence of the correlation function on the SFR 
parameter A. 



VI. SUMMARY AND OUTLOOK 

In this paper we, for the first time, included the long-range N^LO corrections to the 
3NF in the Faddeev calculations. These novel two-pion-exchange, two-pion-one-pion and 
ring interactions supplemented with the one-pion-exchange and contact terms emerging at 
N^LO represent, presently, the most advanced chiral 3NF. We use this force in the triton 
and the doublet neutron-deuteron scattering length Faddeev calculations to fix the two 
low-energy parameters cr, and ce for different sets of regularization parameters, which cut 
off large-momentum or, equivalently, short-range components in the few-nucleon states. 
While the value of the LEG co remains fairly stable, ce features a stronger sensitivity to 
regularization parameters. It will be interesting whether such a behavior is also seen for 
fits involving the complete N'^LO 3NF. We also studied the individual contributions of the 
various topologies to the triton binding energy. The expectation values of the two-pion- 
one-pion and ring terms turn out to be smaller than the ones of the dominant two-pion- 
exchange 3NF for softer values of the regulator. Generally, all expectation values are found 
to be sizable. As expected, we observe a strong sensitivity of the expectation values to the 
regularization parameters. We also looked at the impact of the used 3NF on the two-nucleon 



16 



correlation function in the triton. 

While our work does not yet correspond to a complete N^LO analysis due to the shorter- 
range contributions and relativistic corrections to the 3NF which are not yet available and 
are still missing in our calculations, it does represent a very important step in this direction 
and provides a proof-of-principle that the very complex operator structure of the 3NF at 
N'^LO can be successfully implemented in few-body calculations. In the future, this study 
should be extended to explore effects of the novel terms in the 3NF in few-nucleon scattering. 
This work is in progress. To complete the analysis of the 3NF at N'^LO the inclusion of full 
structure of shorter-range V^iV-cont' '^2n~cont ^s well as V^ ,^ terms should be pursued. The 
numerical implementation of the new terms can be straightforwardly performed using the 
newly developed aPWD scheme, which is successfully tested for the long-range terms in the 
present study. Finally, it should be emphasized that the present work also opens the way 



for applying the novel chiral nuclear forces in many-body calculations, see e.g. |35l437l| for 
some exciting recent developments along these lines based on N^LO 3NFs. 
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